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ABSTRACT 

We use a sample of 53 massive early-type strong gravitational lens galaxies with well-measured 
rcdshifts (ranging from z — 0.06 to 0.36) and stellar velocity dispersions (between 175 and 400 km s -1 ) 
from the Sloan Lens ACS (SLACS) Survey to derive numerous empirical scaling relations. The ratio 
between central stellar velocity dispersion and isothermal lens-model velocity dispersion is nearly unity 
within errors. The SLACS lenses define a fundamental plane (FP) that is consistent with the FP of the 
general population of early-type galaxies. We measure the relationship between strong-lensing mass 
Miens within one- half effective radius (R e /2) and the dimensional mass variable M^ m = G^ 1 a1 2 {Rel2) 
to be log 10 [M lens /10 n M Q ] = (1.03 ± 0.04) log 10 [M dim /10 n Af Q ] + (0.54 ± 0.02) (where a e2 is the 
projected stellar velocity dispersion within R e /2). The near-unity slope indicates that the mass- 
dynamical structure of massive elliptical galaxies is independent of mass, and that the "tilt" of the 
SLACS FP is due entirely to variation in total (luminous plus dark) mass-to-light ratio with mass. 
Our results imply that dynamical masses serve as a good proxies for true masses in massive elliptical 
galaxies. Regarding the SLACS lenses as a homologous population, we find that the average enclosed 
2D mass profile goes as log 10 [M(<i?)/M dim ] = (1.10±0.09) log 10 [i?/i? e ] + (0.85±0.03), consistent with 
an isothermal (flat rotation curve) model when de-projected into 3D. This measurement is inconsistent 
with the slope of the average projected aperture luminosity profile at a confidence level greater than 
99.9%, implying a minimum dark-matter fraction of /dm = 0.38 ± 0.07 within one effective radius. 
We also present an analysis of the angular mass structure of the lens galaxies, which further supports 
the need for dark matter inside one effective radius. 
Subject headings: gravitational lensing — galaxies: elliptical — surveys 



1. INTRODUCTION 

Elliptical galaxies are simple in appearance, but their 
internal structure is resistant to elementary Newtonian 
deduction because the primary luminous tracers of 
their gravity — stars — move not on "cold" circular orbits 
but on "hot" randomized o rbits with a broad distribu- 
tion in phas e space (e.g.. iBertola fc Capacciolil 119751 : 
iBinnevI [1976L H97l llllingworthl 11977ft . Withdetailed 

1 Based on observations made with the NASA/ESA Hubble 
Space Telescope, obtained at the Space Telescope Science Institute, 
which is operated by AURA, Inc., under NASA contract NAS 5- 
26555. These observations are associated with programs #10174, 
#10494, #10587, #10798, and #10886. 

2 Beatrice Watson Parrent Fellow, Institute for Astronomy, Uni- 
versity of Hawai'i, 2680 Woodlawn Dr., Honolulu, HI 96822, USA 
(boltonSif a.hawaii . edu) 

3 Harvard-Smithsonian Center for Astrophysics, 60 Garden St., 
Cambridge, MA 02138, USA (rwaythOcfa.harvard.edu) 

4 Department of Physics, University of California, Santa Bar- 
bara, CA 93101, USA (tt@physics.ucsb.edu) 

5 Sloan Fellow, Packard Fellow 

6 Kapteyn Astronomical Institute, University of Gronin- 
gen, P.O. Box 800, 9700AV Groningen, The Netherlands 
(koopmansOastro . rug . nl) 

7 Institut d'Astrophysique de Paris, UMR7095 CNRS & Univ. 
Pierre et Marie Curie, 98bis Bvd Arago, F-75014 Paris, France 
(gavazzi@iap . f r) 

8 Jet Propulsion Laboratory, California Institute of Technology, 
4800 Oak Grove Drive, M/S 169-327, Pasadena, CA 91109, USA 
(leonidasOjpl . nasa . gov) 

9 Department of Physics and Kavli Institute for Astrophysics and 
Space Research, Massachusetts Institute of Technology, 77 Mas- 
sachusetts Avenue, Cambridge, MA 02139, USA (burles9mit.edu) 

10 Physics Division, Lawrence Berkeley National Laboratory, 
Berkeley, CA 94720-8160, USA (djschlegel91bl.gov) 



spatially resolved observations of the absorption- 
line kinematics of nearby elliptical galaxies, this 
difficulty can be attacked head-on through the use 
of dynamical modeling of the distr ibution func- 
stars in phase sp ace (e.g.. Ivan der Mare! 
Rix fe White! 119921 : IBertin. Saglia. fc Stiaveill 
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2000; iKronawitter et al.l 120001 : 
Gerhard et al.ll2001llCappellari et alj|2006ft . Such efforts 
require a very detailed level of analysis, and rely on 
assumptions of dynamical relaxation and the choice of 
correct dynamical model families. 

Strong gravitational lensing offers an attractive com- 
plement to dynamical modeling, since it measures the 
total mass within the "Einstein radius" defined by the 
lensed images — a scale that is generally comparable to 
the effective radius of the luminous matter distribution 
in the case of galaxy-scale lenses. These lensing mass 
measurements are extremely robust, with few underly- 
ing assumption s and minimal model dependence (e.g., 
lKochaneklfl991h . The main limitation to strong lensing 
as a probe of elliptical galaxy structure has typically been 
the lack of sufficiently large and uniform lens samples. 
This limitation has been overcome in recent years by 
the Sloan Lens ACS Survey (SLACS: IBolton et alj|200a 
Treu et al. 2006; Koopmans et al. 2006; iGavazzi et al.l 
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have begun to deliver statistically significant samples 
of strong lens galaxies by using systematic search 
techniques. 

Previous studies have taken diverse approaches 
to the use of strong-lensing information for the 
measurement of early-type galaxy mass-density struc- 
ture. St rong-lensing aperture mass es have been 
used by IRusin. Kochanek. fe Keetonl (|2003l ) and 
iRusin fc Kochanekl (|2005f ) to provide ensemble con- 
straints upon the mass-density structure of early-type 
galaxies. Lensing mass constraints have been combined 
in self-consistent fashion with stellar dynamics to con- 
strain the mass structure of lenses on a system- by-system 
basis by the Lenses Structure and Dynamics Survey 
(LSD:[K oopmans fc Tredl2^0ll200llTreu"fc Koopmansl 



2002, 2003, 2004) an d the SLACS Survey ( Paper III). 
Kocha nek (2002) and iKochanek et al.l (|2006l ) have used 
the time delays between multiple lensed quasar images, 
together with a given value for the Hubble parameter 
Hq to measure the local mass slope of the lensing 
galaxies. Constraints on the mass profile of lenses based 
upon resolv ed lensed features have been published by 
iCohn et all (l200ll). iMunoz et all (l200l. IWavth et all 



(j2005lK iDve fc Warrenl (l2005lh iBrewer fc Lewis! (|200fif h 



and IDve et all (|2007h . iBolton et aD (|2007l hereafter 



B07) have used lensing aperture masses in combination 
with central stellar velocity dispersions and measured 
effective radii to demonstrate on an empirical scaling 
basis that the mass-density structure of massive elliptical 
lenses from the SLACS Survey is on average independent 
of mass, a result which we establish and explore in more 
detail in the current work. 

In this work, we employ the mass, light, and velocity 
measurements of the full Hubble Space Telescope (HST) 
Advanced Camera for Surveys (ACS) gravitational lens 
sample from the SLACS Survey, published in Paper V, 
to derive numerous empirical scaling relations. The 
analysis of the initial SLACS sample in Paper III 
allowed us to measure an on-average isothermal mass- 
density structure; the current lens sample is large 
enough to investigate trends in structure with mass 
and velocity dispersion. We frame much of our analysis 
and discussion in terms o f the Fundamental Plane 
scaling relation (FP: e.g.. iDiorgovski fc Davisl Il987l: 
Dressier et aT1ll 987: Bender et al. 1992; R.cnzi ni fc Ciottil 
19931: Guzman et all fl9 93: Pah re et al.l 119951 . Il998bl lal: 



Jprgensen et al. 19961: Ivan Dokkum fc Franxl 11996c 



Ciotti et al.lfl99 6: Kels on et aj1ll997l : I Graham fc Collessl 
iScodeggio et al l 119981 : IKochanek et aT " 



van Dokkum & Stanford 
20031: Ivan de Yen et af 



iTreu et all 1200 ll 120021: iBertin et airT2 002 



120031: iBernardi et al l 
120031 : iTruiillo et all 12001 
van der Wei et all l200l Paper II; B07) and the rela- 
tionships between luminosity, "dynamical mass" , and 
strong-lensing aperture mass that illuminate the struc- 



tural explanation for the "tilt" of the FP. In the SLACS 
sample, we have a large number of early-type galaxies 
distributed across the higher-mass end of the FP, with 
uniform and high-quality measurements of the observ- 
ables of redshift, surface brightness, velocity dispersion, 
and effective radius. In addition to these quantities, we 
have a full set of aperture masses — measurements of 
the total mass within the Einstein ring radius. These 
additional mass measurements add another dimension 
of physical constraint to the SLACS sample that is not 
available for other FP galaxy samples, thus offering the 
opportunity to break some of the degeneracy inherent 
in the physical interpretation of FP analyses. 

This paper is arranged as follows. Section [3 gives an 
overview of our measurements and analysis techniques. 
Sections (3H6] present the various scaling relations defined 
by the SLACS sample, including the relative alignment 
and flattening of projected mass and light distributions 
O, the relation between stellar and lensing velocity- 
dispersion measurements O, the FP and "mass plane" 
relations (jJSJ), and various perspectives on the relation- 
ship between dynamical mass, lensing mass, and lumi- 
nosity (|J6]). In <j7] we make a robust determination of 
the ensemble average radial enclosed-mass profile of the 
SLACS sample. Section [8] provides an itemized summary 
of our results, and $9] provides a concluding discussion. 

All computations in this work assume a general rel- 
ativistic Friedmann-Robertson- Walker cosmology with 
matter-density parameter Qm — 0.3, vacuum energy- 
density parameter Q\ = 0.7, and Hubble parameter 
H a = 70kms" 1 Mpc- 1 . 

2. OVERVIEW OF SAMPLE, MEASUREMENTS, AND 
METHODS 

We consider the sample of 63 "grade- A" strong gravita- 
tional lenses presented in Paper V. For analyses employ- 
ing stellar velocity dispersions, we restrict the sample to 
the 53 early-type lens systems with a median SDSS spec- 
tral signal-to- noise ratio of 10 or greater per 69kms _1 
pixel over the rest-frame range of 4100 to 6800A. The 
measurements upon which the current work is based are 
all presented in Paper V, and are described only briefly 
here. 

We use stellar velocity dispersion s measured from 
Sloan Digital Sky Surv e y (SP SS: lYork et all l2000l : 
lAdelman-McCarthv et al.l I2007D spectroscopy, which 
samples a 3"-diameter circular fiber aperture centered 
on the target galaxy. We correct these measured ve- 
locity dispersions (which we den ote by ctsdss) us i ng the 
empirical power-law relation of |j0rgensen et al.l (|1995f ) 
to uniform physical apertures of either R e /8 (to give a 
"central" velocity dispersion cr e g) or R e /2 (o~ e 2, for a spa- 
tial aperture more closely matched to that of the spec- 
troscopy and of the strong-lensing features). We note 
that these corrections are quite small — RMS of 2.5% 
for R e /2 — and that the coefficients of logarithmic scal- 
ing relations we derive are not sensitive to the partic- 
ular choice of velocity-normalization aperture. We use 
magnitudes, effective radii, projected axis ratios, and 
position angles for the lens galaxies measured from fit- 
ting de Vaucouleurs models to ACS-WFC F814W imag- 
ing. To convert from observed magnitudes to rest- 
frame luminosities, we firs t correct for Galactic dust ex- 
tinction using the maps of lSchlegel. Finkbeiner. fc Davisl 
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TABLE 1 

Adopted errors in measured data values 



Measured parameter 


Adopted error 


log 10 Luminosity Ly 


0.010 dex 


log 10 Effective radius R e 


0.015 dex 


log 10 Aperture mass M 


0.010 dex 


iogjg Surface brightness I e = Ly/(2nR^,) 


0.020 dex 


log 10 Velocity dispersion a 


0.030 dex 



(1998), and apply fc-corrections and evolutionary cor- 
rections. Our strong gravitational-lens models provide 
for each lens a robust measurement of the mass Msm 
enclosed within the physical Einstein radius i?Ein (de- 
rived from the angular Einstein radius b together with 
redshifts and cosmology), for both singular isothermal 
ellips o id (S IE: Kas siola fe Kovnerlll993l : iKormann et al.1 
1994; Kc eton fc Kochanekl I1998F and light-traces-mass 
(LTM) models for the lens-galaxy mass-density profile. 
For much of our analysis, we evaluate the strong-lensing 
aperture masses for a uniform physical aperture of i? e /2 
(one-half the effective radius of the luminous compo- 
nent), chosen to be fairly closely matched to the sam- 
ple median Einstein radius i?Ein of approximately 0.6i? e . 
The evaluation of the aperture mass at a radius other 
than i?Ein introduces some dependence upon the as- 
sumed lens-mass model; we gauge this effect by using 
aperture masses from both SIE and LTM models, by 
checking for systematic residual correlations with the 
ratio i?Ein/-Re (of which we find none), and by verify- 
ing that our results do not change when derived using 
only the half of the sample with the smallest fractional 
aperture-mass difference between i?Ein and R e /2. For 
an overview of the theory, phenomenology, and scientific 
applications of str ong gravitational lensing, we refer the 
reader to Part 2 of lSchneider. Kochanek. fc Wambsganssl 
<f2006h . 

Unless otherwise noted, all scaling relations are fitted 
as linear and plane relationships in logarithmic space. 
We define the best fit as that which minimizes the to- 
tal squared orthogonal distance from the line or plane to 
the set of sample data points. Before carrying out the 
fits, we scale all logarithmic data coordinates by their 
typical logarithmic errors, so as to apply a roughly uni- 
form error metric across the multiple observables that 
span the spaces under consideration. The adopted errors 
are taken from the comparisons among multiple measure- 
ment techniques in Paper V, and are given in Table [TJ 
We estimate errors on our fitted paramete rs using boot - 
strap re-sampling of the analysis sample (|Efronlfl979f ). 
Our reported errors are the square-root diagonal entries 
of the parameter covariance matrices constructed from 
the sets of bootstrap parameter estimates. 

3. THE COMPARATIVE ANGULAR STRUCTURE OF MASS 
AND LIGHT 

In addition to the measurement of aperture masses, the 
simplest strong-lens models also give a determination of 
the angular structure of the lens mass profile, in the form 
of either a mass axis ratio (for the SIE models) or an ex- 
ternal shear magnitude (for the LTM models), along with 
the associated position angle. The implications of these 
measurements were discussed at length in Paper III; we 
briefly revisit the analysis here for the larger sample of 



lenses presented in Paper V. 

First we examine the position-angle alignment between 
the major axes of the light profile (as determined by the 
de Vaucouleurs fit) and the mass profile (as determined 
from the SIE model) . For the 58 lenses with light-profile 
axis ratios q s tars < 0.95, the mean position-angle differ- 
ence between mass and light profiles is 2.5° ± 2.4°. The 
RMS alignment is 18°, as compared with an intrinsic er- 
ror of 7.5° in the mass position angles from the analysis 
of Paper V. We can further restrict our attention to a 
subset of 29 of these lenses whose lensed images extend 
over a significant azimuth about the lens-galaxy center, 
for which the mass PA measurement is more tightly con- 
strained, taken from the 32 lenses of the "ring subset" 
described in Paper V. In this case, we find a mean mass- 
light PA difference of 1.4° ±1.9°. From Paper V, the 
RMS error in the mass position angle measurement for 
the ring subset is approximately 2°. The RMS position- 
angle difference for the ring subset of lenses is 10°, equal 
to the value found in Paper III for the original SLACS 
sample. 

Interestingly, the LTM plus external shear models of 
Paper V show a significant amount of preferential align- 
ment of the "external" shear with the major-axis PA 
of the de Vaucouleurs surface-brightness model. Of 
the above-mentioned 29 well-constrained ( "ring subset" ) 
lenses with <7 s t a rs < 0.95, 17 (59%) have their shear PA 
aligned to within ±15° of their luminosity PA, an angu- 
lar range that only encompasses 17% of the meaningful 
range of variation (±90°). If the shear were external and 
uncorrelated with the PA of the luminous distribution, 
this level of random alignment would have a probability 
of less than one in 10 6 . This result suggests that much 
of the shear required by the LTM lens models is not of 
environmental origin, but is instead compensation for an 
intrinsic shortcoming of the LTM lens models. The fact 
that the shear is aligned with the light position angles 
rather than being anti-aligned (90° out of phase) shows 
that the lensed features require a stronger gravitational 
quadrupole at the locations of the lensed images than can 
be provided by the pure LTM models. This translates 
into a requirement either that the radial mass-density 
profile be less centrally concentrated than the light pro- 
file, or that the projected mass profile be flatter than the 
projected luminosity profile. In cither case, the improb- 
able alignment of "external shear" with the major axes 
of LTM lens models provides evidence against the light- 
traces-mass hypothesis and thus in favor of some form of 
dark-matter component. 

We also examine the relative flattening of mass and 
light as measured from SIE and de Vaucouleurs models. 
For the 57 grade-A lenses with early-type morphology, we 
find (gsiE/gstars) = 1-02 ± 0.02 with an RMS deviation 
of 0.12 about the mean: i.e., mass and light have essen- 
tially the same projected axis ratio. As in Paper III, we 
see a trend from mass being rounder than light at lower 
masses towards mass being flatter than light at higher 
masses. Considering the subset of grade-A lenses with 
early-type morphology that both (1) are in the angu- 
larly well-constrained "ring subset", and (2) have mass 
and light position angles aligned to within 15°, this trend 
can be seen as a linear correlation coefficient of qsie / fears 
with (tsie of r — —0.665, corresponding to a significance 
of 99.95%. 
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4. VELOCITY DISPERSIONS: STELLAR AND LENSING 

For SIE lens models, we may translate the measured 
Einstein radii b into lens-model velocity dispersions asm 
through the relation 



1.4 [ 



IE 



47T- 



; SIE 



(1) 



where -Dls and D$ are "angular-diameter distances" 
from lens (foreground galaxy) to source (background 
galaxy) and observer to source respectively. Though the 
conversion is strictly appropriate only for an isotropic, 
spherically symmetric, single-component density distri- 
bution with p oc r~ 2 (in 3D), it provides a robust 
velocity-scale measurement with extremely small statis- 
tical error, to be compared with the (usually noisier) 
central velocity dispersion altars of the stars within the 
lens gal axy. By conside ring bulge-plus-isothcrmal-halo 
models. iKochan ck (1994) suggested that an approximate 
equality should hold between the central stellar line-of- 
sight velocity dispersion and the velocity-dispersion pa- 
rameter of the isothermal halo. This equality has gen- 
erally been confirmed in successively larger samples of 
lens galaxies with mea sured stellar velocity dispersions 
(jKochanek et al.ll2000l LSD, Papers II and III), and the 
current SLACS lens sample allows the most significant 
determination to date. 

Weighting with the statistical errors of the SDSS ve- 
locity dispersions (limited to a minimum error-bar size of 
£cr stars = 0.05a s tars), we find cr SD ss = (0.948 ±0.008)ct S ie 
and <7 e 8 = (1-019 ± 0.008)<7sie- Thus we see that the 
scaling a e s = ogiE holds to within less than 3%, and is 
essentially consistent with being exactly true on average 
at the current level of evidence. The RMS scatter about 
the relation is approximately 0.09 for ctsdss and 0.1 for 
<7 e 8 ) somewhat larger than seen in Paper II. The reduced 
X 2 of this scatter when weighted by the errors as in the 
averaging is 2.7, showing evidence for intrinsic RMS scat- 
ter at the level of about 20kms _1 , or alternatively at 
the level of 7.5% asm (i.e., the necessary value added 
in quadrature to the velocity errors to give a reduced 
X 2 of approximately unity) . Figure [1] shows the ratio 
/ = erstars/csiE as a function of ctsie- Note that our def- 
inition of / follows the conventi on of Paper II, but is th e 
inverse of the definition used in lKochanek et aTl (|2000l ). 

The residuals about the mean cr s tars-versus-(7siE value 
are not significantly correlated with asm, mass, lumi- 
nosity, mass-to-light ratio, effective radius, or ratio of 
Einstein radius to effective radius. There is an apparent 
correlation with cr s tars, but this is an artifact of the rel- 
atively large statistical errors in that quantity. The lack 
of a correlation between / and lens mass in particular 
indicates that the mass-dynamical structure of massive 
early- type galaxies is independent of mass. Likewise, the 
lack of correlation between / and the ratio of Einstein ra- 
dius i?Ein to effective radius R e indicates that the univer- 
sal density structure is nearly isothermal independently 
of dynamical modeling. These two results will be derived 
quantitatively in $6] and [J7] below, and discussed further 
i n m 

Auger (2007) has reported a significant correlation for 
the 15 SLACS lenses analyzed in Paper III between the 
logarithmic mass-profile slopes (from a combined lensing 
and dynamical analysis) and the number of near neighbor 
galaxies to the lenses. The (small) intrinsic scatter of 
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Fig. 1. — Ratio / of stellar velocity dispersion Cstars to velocity- 
dispersion parameter ctsie of the best-fit lens model, as a func- 
tion of 0"siE- Points are shown for altars as observed within the 
3"-diameter SDSS fiber, and as corrected to a uniform physical 
apertu re of R e /8 using the empirical formula of j0rgcnscn et al. 
(1995). For visual clarity, error bars (gray) are only plotted for the 
fiber-aperture points. 

the individual mass slopes about the average value is 
directly related to the scatter in the a e ?, -asm relation ; 
it remains to be determined whether the lAugerl (2007) 
result is confirmed by the larger current SLACS sample 
(Treu et al., in preparation). 

5. THE FUNDAMENTAL PLANE AND MASS PLANE 

The fundamental plane (FP: lDiorgovski fc Davislll987l ; 
iDressler et al.l [l987h is the name given to the approxi- 
mately two-dimensional manifold defined by early-type 
galaxies within the three-dimensional space spanned by 
the logarithms of effective radius R e , surface brightness 
I e , and central velocity dispersion a c . When the FP is 
expressed in the form 

log 10 R e = a log 10 a c + b log I e + d , (2) 

the coefficients a and b are found to be in the approximate 
ranges 1.1 to 1.6 and —0.75 to —0.8 respectively, depend- 
ing upon the sample, wave-band, and me thods of obser- 
vatio n and fitting (see Tables 2 and 4 of iBernardi et al] 
2003). As has been discussed extensively in the litera- 
ture, the implications of this relationship can be under- 
stood in terms of a simple dimensional analysis of the 
form 

(3) 



R„ = cG^alT- 1 !: 



where T is the total mass-to-light ratio within some phys- 
ical aperture and c is a dimensionless constant parame- 
terizing the details of mass-dynamical structure. Equa- 
tion [2] with the values (a, b) — (2,-1) is often referred 
to as the "virial plane" , and the observational fact that 
a < 2 and b > —1 is referred to as the "tilt" of the FP 
relative to the virial plane. The terminology is perhaps 
misleading: the tilt of the FP does not imply a lack of 
virial equilibrium in the individual galaxies that define 
the plane; it simply requires that c and/or T vary as a 
function of the other observables. 

B07 constructed the FP of SLACS lenses and showed 
that the SLACS FP coefficients were consistent with 
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those of the larger SDSS early-type galaxy population, 
while Paper II showed consistency between the SLACS 
lenses and the FP of local galaxies, corrected for lumi- 
nosity evolution. We now revisit the calculation of B07 
for the larger current sample of lenses (53 early-type 
lens galaxies with well-measured SDSS velocity disper- 
sions), using the data presented in Paper V. 11 We de- 
scribe the FP in the form of Equation [2 with <j e 2 tak- 
ing the role of central velocity dispersion a c . As noted 
above, the aperture of R e /2 is chosen to reflect the ap- 
proximate median ratio of Einstein radius i?Ein to R e 
within the sample. Due to the velocity- dispersion aper- 
ture correction formula that we employ (jjOrgcnsc n et all 
|l995[h the er e2 values are related trivially to the more 
commonly quoted a e g, values by a fixed factor, with 
logio cr e8 = lo Sw a e2 + 0.0241. The best-fit FP coeffi- 
cients as determined by the method described in are 
presented in Table [2] The plane is shown in edge-on 
projection in the left panel of Figure [2] The RMS resid- 
ual scatter in log 10 R e of the fit is 0.064 dex. The RMS 
error-scaled orthogonal residual scatter of the fit is 1.44, 
reasonably consistent with error estimates but indicat- 
ing a degree of intrinsic scatter. The bes t-fit coefficients 
are co nsistent with those determined by iBernardi et all 
(2003) from SDSS data for a larger early-type galaxy 
sample, though the details of that analysis are somewhat 
different than ours. 

As originally noted in B07, there is a significant corre- 
lation between log-radius residuals about the best-fitting 
FP and the rest-frame F-band mass-to-light ratio Ty 
(for mass and light within R e /2, evaluated from B-spline 
light models and SIE mass models). This correlation has 
a linear coefficient r = —0.337 (corresponding to a 98.6% 
significance) and is in the sense that lens galaxies with 
effective radii larger than their FP-predicted values tend 
to have lower mass-to-light ratios. This is to be expected 
from a simple dimensional analysis as in Equation [3J if 
er and / are held fixed (i.e., for a given point on the FP), 
an increase in R should give a decrease in Ty. This cor- 
relation suggests observationally that the scatter about 
the FP is due at least in part to an intrinsic scatter in 
mass-to-light ratios at a given point in the plane. 

In addition to the FP, we consider the "mass plane" 
(MP) of SLACS lenses as defined in B07. By using 
strong-lensing aperture masses (corrected to the uniform 
aperture of i? e /2), we can replace the surface brightness 
I e with the surface mass density £ e 2 within R e /2. The 
MP formalism is attractive in that it is independent of 
any luminosity-evolution effects that occur in a spatially 
uniform manner within the lens galaxies. Thus the MP 
can be tracked across redshift without correction for the 
dimming of stellar luminosity. The MP should also be 
better suited to comparison with theory and numerical 
simulation of galaxy formation, merging, and dynamical 
evolution processes. Since the scale length of the MP is 
still taken from the luminosity distribution and the MP 
velocities are still traced by the stars, it is not a pure mass 
space formulation. However, one can extract total (lumi- 
nous plus dark matter) aperture masses, stellar effective 
radii, and stellar velocity dispersions from any sensibly 

1 Our current determination also differs from that of B07 in the 
use of luminosities based upon de-bugged fc-corrections as described 
in Paper V. The overall conclusions of B07 are unchanged. 



TABLE 2 
FP AND MP COEFFICIENTS 



Surface 
Density 


a 


b 


d 


Ie 


1.28 ± 0.22 


-0.77 ±0.07 


-0.09 ± 0.07 


£e2,SIE 


1.82 ± 0.19 


-1.20 ±0.12 


0.91 ±0.14 


S e 2,LTM 


2.10 ± 0.24 


-0.83 ±0.12 


0.56 ±0.15 



Note. — Coefficients are for FP and MP relations in the form 
log 10 R e = a log 10 <r e 2 + blogjQ (surface density) + d. R e measured 
in kpc, a e 2 in units of 100 kms -1 , I e in units of 10 9 Ly q/ kpc - 2 , 
and S e 2 in units of IO^A/q/ kpc -2 . Units are chosen so as to 
reduce artificial covariance between plane coefficients and plane 
zero-points. Fit method is described in J2] 



constructed theory of galaxy formation, evolution, and 
merging. One cannot extract luminosities without delv- 
ing into the entirely different domain of stellar popula- 
tions and their evolution. Thus the MP formulation — 
and more elementally, the data upon which it is based — 
represent an important step in bringing theory and ob- 
servation together by separating the subject of stellar 
populations from the subject of galaxy formation. 
We express the MP in a form analogous to the FP: 



log 10 R e = a m log 10 cr e2 + b m log 10 S e2 + d 



(4) 

Since E e 2 depends somewhat upon the lens model used 
to evaluate the surface density, we fit for MP coefficients 
using E e 2 values from both the SIE (singular isothermal 
ellipsoid) and LTM (light-traces-mass) mass models of 
Paper V. Values for the best-fit coefficients are given 
in Table H The best-fit MP (using the SIE aperture 
masses) is seen in edge-on projection in the right-hand 
panel of Figure O The RMS residual scatter in log 10 R e 
about the best-fit MP is 0.076 dex using SIE surface den- 
sities and 0.10 dex using LTM surface densities. The 
RMS error-scaled orthogonal residuals are 1.24 (SIE) and 
1.50 (LTM). For the case of the SIE-based MP, the plane 
is tighter than the FP when judged by the RMS error- 
scaled orthogonal residuals (the quantity minimized in 
the fitting), but not exceptionally so: the edge-on MP is 
not noticeably tighter than the edge-on FP as seen in of 
Figure [2j 

We can see from Table [2] that the coefficients of the 
MP are in a sense less "tilted" relative to the values 
(a m ,b m ) — (2,-1). To quantify this, we consider 10 6 
bootstrap re-samplings of our lens sample, from which we 
compute confidence limits in the a m -b m plane. The val- 
ues (a m , b m ) = (2, —1) fall on the contour enclosing 77% 
(60%) of the re-sample coefficient points for SIE (LTM) 
based fits, showing that the MP is in fact fairly consis- 
tent with the plane represented by {a m ,b m ) = (2,-1). 
We cast this statement in more physical terms in Sj6j 

6. LENSING MASS, DYNAMICAL MASS, AND LUMINOSITY 

Depending on the details of the coefficients, the FP 
can be understood in terms of a systematic relation be- 
tween luminosity and a dimensional mass variable 12 such 
as M dim = G-V 2 9 (j? e /2) of the form L oc M," : (e.g 



iFaber et al.lll987t iBender et alll992t ICiotti et al.lll996h . 
Such a relation leads to FP coefficients in Equation [5] 

12 Our "dimensional mass" Af(ji m is, modulo a constant factor, 
equal to the "effective mass" or "dynamical mass" referred to by 
other authors. We use the term "dimensional" to avoid confusion. 
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Fig. 2. — Edge-on projections of the best-fit FP (left) and MP (right) relationships for the SLACS lens sample. Units are 100 km s 1 
for (T e 2, 10 9 Lyg kpc — 2 for I e , lO 9 M0kpc -2 for S e 2, and kpc for R e . These are the same units used in the FP and MP fits, and were 
chosen so as to lessen covariance between plane coefficients and plane zero-point. Diagonal solid lines represent equality between log 10 R e 
as observed and as predicted from the best-fit FP. 



given by a = 2rj/ (2 — if) and b = —1/(2 — 77). This may 
be conceptualized as a systematic variation of a "dimen- 
sional mass-to-light ratio" T<ji m = Mdim/L with Mdi m . 
It is, however, important to keep in mind that Ydi m is 
linearly proportional to the true mass-to-light ratio T 
only if the FP tilt is due to a systematic mass-to-light 
ratio trend, rather than to a trend in mass-dynamical 
structure. 

The availability of strong-lensing aperture masses in 
addition to the traditional FP observables allows us to 
directly test the alternative hypotheses for the "tilt" of 
the FP, because in addition to M<ji m , we measure the 
aperture mass Mi ens from strong lensing. Consider the 
following two relations: 

Lv = ciM^ m , (5) 
M lens = coM^ im . (6) 

Consider also the following two alternative hypotheses: 
(1) on average, early-type galaxies have a universal mass- 
dynamical structure, and the tilt of the FP is due to a 
systematic trend in total mass-to-light ratio with mass; 
or (2) on average, early-type galaxies have a universal 
mass-to-light ratio, and the tilt of the FP is due to a sys- 
tematic trend in mass-dynamical structure. If hypoth- 
esis (1) is correct, then we should find 6 ~ 1, whereas 
if hypothesis (2) is correct, we should find 6 ~ 77. Of 
course, we may also find that r\ < 8 < 1, since in princi- 
ple both explanations could contribute to the tilt of the 
FP. Roughly speaking, hypothesis (1) represents "ho- 
mology" and hypothesis (2) represents "non-homology" . 

Before proceeding, we note that numerous works 
have explored the possible role of a systematically 
varyi ng Sersic index n in causing the tilt of the FP 
(e.g.. iHiorth fc Madsenl fl995l iCraham fc Collessl Il997t 
iBertin et al.ll2002tlTruiillo et al.ll2004h . and thus we must 
not disregard this possibility in our own analysis. We 
compute n for all SLACS lenses by continuing the de 
Vaucouleurs model optimizations after freeing the index 
from its fixed n = 1/4 value. We find that n is com- 



pletely uncorrelated with lensing mass, dynamical mass, 
luminosity, and velocity dispersion within the sample, 
and therefore the inclusion of the Sersic index as a sig- 
nificant factor in our analysis is not motivated by the 
data. This lack of correlation between n and other quan- 
tities is in fact consistent with other studies, since the 
SLACS sample is confined to relatively high-mass/high- 
luminosity early- type galaxies, and does not extend over 
a sufficient range to define the n-L correla tion signifi- 
cantly given the level of intrinsic scatter (e.g..lCaon et al.l 
1991 iD'Onofrio et all H99l iGraham fc Guzman! 120031 
Ferrarese et al.ll2006f ). 

Table [3] shows the results of fits for the normalizations 
and exponents of the relations defined in Equations O 
and[6j as well as for the relation between luminosity and 
lensing mass. Within the uncertainties, the clear result 
is that S ~ 1 while r\ < 1: thus, we conclude based upon 
our sample of lenses that the tilt of the FP — as defined by 
massive ellipticals — is due to a systematic trend in mass- 
to-light ratio and not to a trend in mass-dynamical struc- 
ture. These relations are shown graphically in Figure [31 
where we can see by eye that the logarithmic slope of the 
-Miens-versus-Mdim relation is significantly less shallow 
than those of the Ly-versus-Mdi m and Ly-versus-Mi ens 
relations. The slopes of the L-versus-M relations change 
only negligibly when the de Vaucouleurs model lumi- 
nosity is replaced with the more general radial B-spline 
model aperture luminosity within R e /2 (described in Pa- 
per V). If we fix 5 = 1, we find (log 10 [Mi cns /M dim ]) = 
0.530 ±0.012 (RMS of 0.08) for SIE aperture corrections 
and (log 10 [Mie„ s /M dim ]) = 0.543 ± 0.014 (RMS of 0.1) 
for LTM corrections. Accounting for our best estimates 
of the typical measurement errors, the intrinsic scatter 
in log 10 [Mispa/Msm] about the mean value that gives 
a reduced % 2 of approximately unity is 0.057 dex RMS, 
or about a ±13%. This is consistent with the intrinsic 
velocity-dispersion scatter of 7.5% seen in given that 
Mdim scales as the square of velocity dispersion. 

The result that S ~ 1 while -q < 1 is one of the central 
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TABLE 3 
Mass and light scaling relations 



Scaling Relation Prefactor Exponent 

(Msie/IC^A/q) = coCMdim/lO 1 ^-/©) 4 log 10 c = 0.54 ± 0.02 S = 1.03 ± 0.04 

(A^ltm/1O u M ) = c o (M dim /lO 11 M )' 5 log 10 c = 0.54 ± 0.02 S = 0.99 ± 0.05 

(L v /l0 n L a ) = ci(M dim /lO 11 M )'' log 10 ci = 0.16 ±0.02 r\ = 0.77 ± 0.04 

(L v /W n L @ ) = c 2 (M SIE /10 11 M Q ) r '' log 10 c 2 = -0.24 ± 0.01 r/ = 0.73 ± 0.03 

(L V /1O 11 L ) = c 2 (M LT m/1O 11 M )'' / log 10 c 2 = -0.26 ± 0.01 rj = 0.77 ± 0.03 
Note. — Mass values AfsiE an d JWlxm are niasses within R e /2 as evaluated from SIE and LTM lens models. Dimensional mass variable 
is defined as Af dim = G _1 o-^ 2 (^?e/2). Fit method is described in fj2] 




Fig. 3. — Relationships between mass and luminosity variables. Left panel: Lensing mass Afi ens within one-half effective radius (evaluated 
from the SIE lens models) versus dimensional mass variable Af d ; m = G~ 1 <j^ 2 (Re / 2) ■ Center panel: Rest- frame V-band luminosity L versus 
A^dim- Right panel: Rest-frame V-band luminosity L versus Afi cns . In each panel, the best-fit linear relationship between the two 
logarithmic quantities is shown by the solid line. Dashed lines indicate a slope of unity, for reference. Typical errors in Af d ; m are shown by 
gray error bars in the left and center panels; errors in Afi ens and L are smaller than the plot symbols. 
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Fig. 4. — Distribution of n and <5 parameter fits to 10 6 bootstrap 
re-samples. Results are shown for lensing aperture masses eval- 
uated using both SIE lens models (black) and LTM lens models 
(gray). Dashed lines of S = 1 and of 5 = r\ are shown as well. 

results of this work and of B07, and thus its significance 
merits special attention. Specifically, we consider the 
distribution of pairs of r\ and 5 values fitted to the same 
sets of bootstrap samples. For both the SIE and LTM 
aperture masses, we generate 10 6 bootstrap re-samples 
and fit rj and 5 for each. The resulting distributions can 
be seen in Figured] We see that for both mass models, 5 
is significantly greater that r/ and consistent with 5=1, 
while rj = 1 and 6 = rj are ruled out at high significance. 



The skeptic might worry that our result of a universal 
mass-density structure as a function of mass has been 
"baked in the cake" by our use of a single universal pro- 
file (i.e., the isothermal model or the approximately dc 
Vaucouleurs LTM model) to evaluate aperture masses 
for the uniform physical aperture of R e /2. The robustly 
measured quantity is the m ass within the Einstein ra- 
dius i?Ein (|Kochanekj 1 1 99 lh . and in general i?Ein differs 
from R e /2 (though, by design, not by a large amount). If 
such an explanation were masking a true mass-dependent 
trend in structure, we would expect to see a correlation 
of REin/Re with the residuals about the best-fit Mi ons - 
versus-Mdi m relationship. In fact, we see no such correla- 
tion. In addition, if we restrict our analysis to that half of 
the lens sample for which i?Ei n is most nearly matched to 
R e /2 (as quantified by the fractional difference between 
the aperture masses within i?Ein and R e /2), we see no 
significant changes in our results. 

We can easily translate our relations between V-band 
luminosity and lensing mass within R e /2 (the last two 
rows of Table [3]) into expressions for the mass-to-light 
ratio within R e /2 as a function of luminosity. Using the 
fact that 32.0% of the total de Vaucouleurs model flux is 
enclosed within the R e /2 aperture, we have 



V,e2 



-1/77/ 



1 



0.320 



u 



(7) 



For SIE aperture masses, we find log-^Ty^/Tv^) = 
(0.83 ± 0.01) + (0.37 ± 0.06) log^i^/lO 11 ^©), while 
for LTM aperture masses we obtain log 10 (Tv je 2/Ty j Q) = 
(0.83 ± 0.01) + (0.29 ± 0.05) log lo (Lv/lO n Lv, )— i.e., a 
total central V-band mass-to-light ratio of approximately 
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6.8 times solar for a V-band luminosity of 10 11 times so- 
lar. This T-L relation should not necessarily be regarded 
as fundamental, however: we find comparably significant 
correlations of T with Mi ens , Mdi m , R e , and a e %. 

The details of our mass-luminosity relations are sensi- 
tive to the assumed rate of passive evolution which we 
have attempted to remove from the sample. The SLACS 
lens sample exhibits a significant luminosity-rcdshift de- 
generacy (see Paper V) which prevents us from solving 
simultaneously for the mass and redshift dependence of 
luminosity The main effect that we expect is a cer- 
tain level of unmodeled luminosity dependence in the 
rate of V^-band evolution, characteristic of more recent 
star formation (and hence faster fading) in less lumino us 
galaxies (e.g.. ICowie et al.l Il996t ITreu et al.l l2005allbl ; 
Ivan der Wei et al.ll2006f ). though the effect should be less 
pronounced in the V band than in the B band. In the ab- 
sence of published constraints on the mass dependence of 
y-band luminosity evolution, we consider the effects of a 
simple "toy model" whereby the rate of luminosity evolu- 
tion varies from d log 10 Lyjdz = 0.6 to 0.2 linearly over 
the range log 10 (M lens /M Q ) = 10 to 12. (This is likely 
to be an e xtreme scenario — -s ee the rest-frame -B-band 
analysis of ITreu et aT1l2005bl ) In this case, the expo- 
nents of the L- versus- M relations given in Table [3] change 
from j) = 0.77 to X] = 0.82 and rf = 0.73 to rj' = 0.78, 
with unchanged uncertainties. We defer more detailed 
analysis until the completion of multi-band observations 
currently underway with _H5T-WFPC2, which will en- 
able the determination of rest-frame £?-band luminosi- 
ties and subsequent connection to a mo re comprehensive 
literature on luminosity evolut i on fe.g.lTreu et al.ll2002l 
2005bt Ivan der Wei et aTll2005l : Idi Serego Alighieri et all 
2005). These ambiguities of luminosity evolution have 
no bearing upon the pure mass-dynamical relations, and 
no effect upon our result that the structure of massive 
elliptical galaxies is (on average) independent of galaxy 
mass. 

7. ROBUST ENSEMBLE MEASUREMENT OF THE RADIAL 
MASS-DENSITY PROFILE 

In this section, we constrain the average projected 
mass-density profile of our lens galaxies without any dy- 
namical modeling by considering the lens ensemble as 
a whole, and using the fact that the aperture masses 
are measured most robustly for physical apertures (i.e., 
Einstein radii) ranging from 0.2 to 1 times the effec- 
tive radius depending upon the individual lens. 1 Mo- 
tivated by the FP, we regard the sample as a two- 
parameter family, parameterized by effective radius R e 
and velocity dispersion o~ e 2- From these quantities, 
we construct the previously defined "dimensional mass" 
Afdim = G _1 cTg 2 (i? e /2). We then scale the measured 
aperture masses within the Einstein radii M(< R) — 
-^Ein by .Mdim, and scale the corresponding apertures 
R = i?Ein by the effective radii R e . The resulting di- 
mensionless projected mass-radius relation is shown in 
logarithmic coordinates in Figure [5] Describing this re- 
lation with a power law (i.e., linear in the logarithmic 

13 In fact, the sample probes three-dimensional radii even be- 
yond the physical scale of R e , due to the sensitivity of lensing to 
projected mass. For the singular isothermal sphere, 36% of the 
mass (the exact fraction is 1 — 2 /it) within a cylinder of radius R 
is exterior to the corresponding sphere of the same radius. 



coordinates) in the form 

log 10 [A/(<i?)/A/ dim ] = g log 10 [R/R e ] +h, (8) 

we find g = 1.10 ± 0.09 and h = 0.85 ± 0.03 for SIE aper- 
ture masses or g = 1.13 ± 0.11 and h = 0.85 ± 0.03 
for LTM aperture masses. RMS orthogonal residuals 
scaled by error estimates are 1.27 (SIE) and 1.41 (LTM), 
indicating reasonable consistency of the data with the 
fitted relation, but with evidence for intrinsic scatter 
as in other relations. The residuals about the best- 
fit relation are completely uncorrelated with Mdi m , an 
echo of the mass- independent structure result of Sj6] For 
a pure singular isothermal model, the projected aper- 
ture mass scales as M(<R) cx R — i.e., g = 1. Thus 
we find that our data are consistent with the univer- 
sal total-mass profile being isothermal. Translating back 
into a three-dimensional power-law profile, our results 
correspond to p(r) oc r - 1 - 90 ± 09 (using SIE aperture 
masses) and p{r) oc r ^ 1 - 87 ± 010 (using LTM aperture 
m asses). This a nalys is is very similar to that emplo yed 
by IRusin et"all (|2003h and lRusin fc Kochanekl (|2005D to 
derive constraints on the quasar-lens population. The 
key difference for our application to SLACS is that, since 
we have stellar velocity dispersions for our entire sam- 
ple, we can directly scale lensing masses by dynamical 
masses, ignoring issues of luminosity trends and evolu- 
tion rates . The most likely averag e power-law profile 
found by iRusin fc Kochanekl {2005) is slightly steeper 
(p(r) oc r - 2 06 ± 017 ) than that found in this work, but 
the two results are consistent within their combined er- 
rors. Our (nearly) isothermal result is also consistent 
with the lens+dynamical modeling resul ts of LSD and 
Pape r III, with the dynamical analysis of iGerhard et al.l 
120011 and with the strong+weak lensing analysis of Pa- 
per IV. The particular advantage of the method pre- 
sented here is that it is simple and robust, with little 
susceptibility to assumptions or system atic errors, (e.g., 
IGerhard et aIH200lt IRusin et al.ll200l LSD; Paper III; 
Paper IV). 

For comparison, the logarithmic slope of the projected 
enclosed luminosity profile changes very little over the 
range 0.3 to 1 times R e where the bulk of our data are 
concentrated, having a value of gu g ht — 0.71 for the de 
Vaucouleurs model and a mean value of <?ii g ht = 0.68 
for the more general B-spline luminosity profile models 
of Paper V. If a light-traces-mass model were appro- 
priate, we would expect to find similar values for the 
logarithmic slope of enclosed mass with radius. The one- 
sided 99%, 99.9%, and 99.99% lower limits on the mass 
slope variable g from bootstrap re-samples are 0.85, 0.76, 
and 0.67 respectively for SIE fits (0.83, 0.72, and 0.62 
for LTM fits). Thus with only three basic observables — 
velocity dispersion, effective radius, and Einstein radius 
(as translated into physical units and enclosed mass using 
spectroscopic redshifts) — and no dynamical modeling, we 
can falsify the light-traces-mass hypothesis at 99.9% to 
99.99% confidence even inside the effective radius, where 
dark matter is often assumed to play little role. This re- 
sult argues for a preference of SIE lens model parameters 
over those from the LTM models. Figure [5] shows the 
slope of the de Vaucouleurs aperture-luminosity profile 
for comparison with the enclosed-mass data points; the 
inconsistency of the data with an LTM model is visually 
apparent. 
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Fig. 5. — Dimensionless projected mass-radius relation. Hori- 
zontal axis is mass-aperture (Einstein) radius scaled by lens-galaxy 
effective radius. Vertical axis is lensing aperture mass within the 
Einstein radius scaled by the dimensional mass variable formed 
from the combination of central velocity dispersion and effective 
radius. Error bars show the typical error in each coordinate used 
in the scaled orthogonal fit for the best relation. The solid line 
shows this best-fit linear relation in the logarithmic coordinates, 
given by log 10 [M(< i?)/M dim ] = (1.10±0.09) log 10 [R/R e ] + (0.85± 
0.03). For comparison, dashed lines show the logarithmic enclosed- 
luminosity profile for the de Vaucouleurs surface-brightness model, 
with various offsets. If a light-traces-mass model were correct, the 
data points would trace the slope of these de Vaucouleurs curves. 

We can translate our measurement into a lower limit 
on the average central dark-matter fraction by making 
a "maximal bulge" assumption — specifically, assuming 
that the total mass and the stellar mass are equal within 
0.3i? e . We then take the projected stellar mass-profile 
to go as M s tars(< R) oc R 07 (following the luminos- 
ity profile) and the projected total-mass profile to go 
as Mtotai(< R) i? 11±01 . This gives a lower limit 
on the average projected dark-matter fraction of /dm = 
0.38 ± 0.07 within one effective radius R e . This is con- 
sistent with the highest values found thro ugh dynamical 
analy sis of nearby elliptical galaxies (e.g., iGer hard et all 
12001 . 

8. SUMMARY OF RESULTS 

By combining the direct mass measurements provided 
by the strong gravitational lensing effect with traditional 
measures of galaxy luminosity, size, and velocity disper- 
sion, we have derived empirical scaling laws for massive 
early- type galaxies. The results of our analyses are sum- 
marized as follows. We remind the reader that these re- 
sults apply specifically to galaxies with velocity disper- 
sions between approximately 175 and 400 km s , rest- 
frame F-band luminosities between 10 10 and 3 x lO 11 !/©, 
and redshifts between 0.06 and 0.36. (The few higher 
redshift SLACS lenses presented in Paper V have an 
SDSS spectral SNR that is too low for a reliable velocity- 
dispersion measurement.) 

1. The projected major axes of the mass and light 
distributions within the SLACS lens sample are 
aligned to (APA) = 2.5° ± 2.4° (RMS of 18°) for 
55 lenses with light axis ratios fears < 0.95. Re- 
stricted to the subset of 29 lenses whose angular 



structure is tightly constrained by lensed images 
through an extended azimuthal range about the 
lens center, this alignment tightens to (APA) = 
1.4° ± 1.9°, with an RMS of 10°. 

2. For light-traces-mass (LTM) lens models, the re- 
quired external shears are preferentially aligned 
with the major axes of the lens galaxies at ex- 
tremely high significance when the sample as a 
whole is considered. This alignment suggests a fal- 
sification of the LTM hypothesis. 

3. The average relative flattening between mass and 
light for the SLACS sample, measured by the ratio 
of projected (2D) minor-to-major axis ratios gsiE 
(mass) and g stars (light), is (gsra/fears) = 1.02 ± 
0.02 (RMS of 0.12): i.e., consistent with unity on 
average. We see a decreasing trend in <?siE/9stars 
with increasing velocity dispersion asm (measured 
from the SIE lens models). 

4. The ratio / = Ces/csiE of central stellar velocity 
dispersions to isothermal lens-model velocity dis- 
persions is / = 1.019 ± 0.008. The RMS scatter 
about the mean is approximately 0.1, which when 
considered along with measurement errors corre- 
sponds to an intrinsic velocity scatter of about 
20kms~ 1 (or about 7.5% of ctsie)- The residual 
scatter about the mean relation is not correlated 
with mass, luminosity, velocity dispersion, mass- 
to-light ratio, effective radius, or ratio of Einstein 
radius to effective radius. 

5. The SLACS lens galaxies define a fundamental 
plane (FP) that is consistent with the FP measured 
for the general population of early- type galaxies ob- 
served by the SDSS and other surveys. 

6. The SLACS lens sample also defines a "mass 
plane" (MP) relation, obtained by replacing sur- 
face brightness with surface mass density as mea- 
sured from strong lensing. The MP is significantly 
less "tilted" than the FP, and is essentially consis- 
tent with simple expectations based on the virial 
theorem and a universal mass-dynamical structure 
within the population. 

7. The relationship between rest-frame V- 
band luminosity Ly and the dimensional 
mass variable Mdim = G~ 1 al 2 (R e /2) is 
given by log lo [iy/lO u i ] = (0.77 ± 
0.04)log 10 [M dim /10 n M Q ] + (0.16 ± 0.02). 
A similar relation holds when Mdim is re- 
placed with the strong lensing-determined 
mass within i? e /2: log 10 [Ly/10 11 L©] = 
(0.73±0.03)log 10 [M lens , e2 /10 11 M Q ]-(0.24±0.01). 
The consistency between the slopes of these two 
relations implies that the systematic "dynamical" 
mass-to-light variation inferred from the FP is 
representative of a true mass-to-light variation 
within the early-type galaxy population. 

8. The relationship between Mdi m and the mass 
M^ns within R e /2 as measured by strong lens- 
ing is given by log lo [M lens /lO n M ] = (1.03 ± 
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0.04) log lo [M dim /lO 11 M ] + (0.54 ±0.02). The uni- 
tary (within errors) slope of this M-versus-M re- 
lation indicates that the mass-dynamical structure 
of early-type galaxies does not vary systematically 
with mass over the range probed by the SLACS 
sample (approximately 9.8 < log 10 [Af d i m /M Q ] < 
11.6). Our result indicates that the dimensional 
mass is a suitable proxy for the true mass within 
the central regions of massive early-type galax- 
ies. Fitting directly for an overall scaling be- 
tween lensing and dimensional masses, we find that 
(log 10 [M lens /M dim ]) = 0.530 ± 0.012 with an RMS 
scatter of 0.08 dex about the mean. After account- 
ing for measurement errors, this scatter is approx- 
imately 0.057dex (or about ±13%). 

9. The difference between the slopes of the M-L re- 
lations of result [7] and the M-M relation of re- 
sult [5] is very highly significant, and implies that 
the "tilt" of the FP is due to a systematic varia- 
tion in the total mass-to-light ratio with mass or 
luminosity, rather than to a systematic variation 
in mass-dynamical structure. The SLACS sample 
shows no significant correlation between luminosity 
and Sersic index n. 

10. Expressed as a function of luminosity, the central 
V-band total mass-to-light ratio Ty, e 2 of SLACS 
lenses scales according to log 10 (Tv >e 2/Ty : Q = 
(0.83 ± 0.01) + (0.37 ± 0.06) log^Ly/lO 11 ^^). 

11. We obtain a nearly model-independent ensemble 
constraint on the average mass-density profile in 
the SLACS sample by assuming (in accordance 
with the FP) that the lenses form a two-parameter 
family. We scale strong-lensing aperture masses by 
Mdi m and radial mass apertures (Einstein radii) 
by i? e , and find a non-dimensional projected mass- 
radius relation given by log 10 [M(< -R)/M dim ] = 
(1.10±0.09)log 10 [i?/i? e ]±(0.85±0.03). This result 
is inconsistent with the steeper slope of the pro- 
jected aperture luminosity profile at a level of sig- 
nificance between 99.9% and 99.99%, thus falsify- 
ing the light-traces-mass hypothesis. If we assume 
that all projected mass interior to 0.3i? e is in the 
form of stars — a "maximal bulge" assumption — 
this result translates into a lower limit on the 
average projected dark-matter fraction of /dm = 
0.38 ± 0.07 inside one effective radius. The three- 
dimensional mass-density profile corresponding to 
our two-dimensional result is p(r) oc r - 1 - 90 ± 09 ) 
consistent with the isothermal (flat rotation curve) 
model. 

9. DISCUSSION AND CONCLUSIONS 

The SLACS lens sample provides a unique resource 
for the quantitative study of the mass-dynamical struc- 
ture of massive early-type galaxies. This is due to the 
addition of strong-lensing aperture masses to a full com- 
plement of traditional galaxy observables over a signifi- 
cant range of intrinsic parameter variation. The fact that 
the SLACS gravitational lens sample defines a fundamen- 
tal plane similar to that defined by the larger sample of 
SDSS early-type galaxies suggests that deductions based 
upon the SLACS lenses can be generalized to early-type 



galaxies in general. This conclusion is further supported 
by the FP analysis of Paper II, and by the analysis in Pa- 
per V of the distribution of SLACS lenses in luminosity 
within their parent samples. 

9.1. The mass-velocity connection 

The lack of correlation between / = CTcs/csie and ei- 
ther mass or R Em / R e in fact contains the essence of our 
results on the mass independence of galaxy structure and 
the near isothermal nature of the radial profile, as we now 
illustrate. Consider the following relation for the mass 
enclosed within the Einstein radius, which is a conse- 
quence of len sing geometry and holds for al l mass density 
models (e.g.. iNaravan fc Bartelmairn]|1996l ): 

M(R Ein ) = M(bD L ) = • (9) 

Furthermore, take Equation [T] as the definition of the 
lens-model velocity dispersion parameter ctsie in terms of 
the observable angular Einstein radius b (which is in turn 
related to the physical Einstein radius through i?Ein = 
bD E ). Now consider an idealized case where i?Ein = -R e /2 
across a range of masses. Equation [9] becomes 

, 2 \ ( Miim\ 

= *"(<*3Ie) T2 

V "stars / 

= vr/- 2 A/ dim . (10) 

Thus the lack of a correlation of / with M& lTCl implies a 
linear relationship between M(R e /2) and Mdim- Simi- 
larly, if we now consider R = i?Ein for any fraction of R e 
and divide Equation [5] by M d i m , we obtain 

M(R) _ fc 2 b D S \(R\( G \ 
M dim K \Att D LS J \Gj W stais R e /2j 
= 2na 2 slE [R/(a 2 stal . s R e )] 

= 2nf- 2 (R/R e ) . (11) 

Thus a non-isothermal mass profile, in which the enclosed 
mass does not grow linearly with radius, would appear as 
a correlation of / with the aperture-radius ratio R/ R e . 

9.2. The FP-structure connection 

The FP is often contrasted with the so-called "virial 
expectation" , but there is in fact no a priori reason to ex- 
pect that either the stellar mass-to-light ratio or the cen- 
tral dark-matter fraction should be constant with galaxy 
mass. Nevertheless, the existence of the FP of elliptical 
galaxies implies a regularity in their formation and evolu- 
tion history. Our result explains this regularity in terms 
of a universal mass-dynamical structure which is the end 
state of massive elliptical evolution regardless of mass, 
together with a systematic trend in total (luminous plus 
dark) mass-to-light ratio with galaxy mass. As discussed 
in $61 there is no apparent trend of Sersic index n with 
luminosity within the SLACS sample. Despite this fact, 
the SLACS lenses define a clear FP relationship; the ex- 
planation of the SLACS FP must therefore lie with other 
factors. We reiterate here, though, that the SLACS lens 
sample is confined to the high-mass end of the elliptical 
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galaxy population. Thus we cannot rule out the impor- 
tance of the Sersic index to the FP of lower-mass early- 
type galaxies. An important implication of our result 
is that the "dynamical mass" can be used as a suitable 
proxy for the true mass inside 1-10 kpc, with a conver- 
sion factor that is independent of gal axy mass or size. 
A sim ilar result has been reported by ICappellari et al.1 
(2006) for their dynamical analysis of a mostly lower- 
mass galaxy sample. 

In the strictest sense, the fact that we find no mass- 
dependent trend in the ratio of dynamical mass to true 
mass is suggestive of a universal mass-dynamical struc- 
ture, but not fully conclusive. In principle, a com- 
bination of various other dynamical details — perhaps 
trends in the anisotropy profile of the stellar orbits, 
perhaps trends in the extent of dynamical relaxation — 
could conspire to give an unchanging dimcnsionless struc- 
ture constant with mass. We also note that the mass- 
redshift degeneracy in the SLACS sample could in prin- 
ciple also be masking a mass-dependent trend if signifi- 
cant dynamical evolution occurs with redshift (see, e.g., 
Ivan der Marel fc van Dokkuml [2007aU H). and the more 
massive lens galaxies evolve to have the same dynami- 
cal structure as the less massive lens galaxies at lower 
redshift. 

9.3. Stellar mass or dark mass? 

Gravitational lcnsing in both strong and weak forms 
measures total mass: stars and dark matter together. 14 
A similar statement holds for mass constraints based 
upon dynamical measurements (which we note here even 
though our use of stellar dynamics in this paper is lim- 
ited to empiricism and dimensional analysis): although 
the dynamical tracers are distributed with the optical 
luminosity, their orbits are determined by the potential 
of the total (luminous plus dark) matter distribution. 
The results presented here strongly suggest that the total 
mass-density structure of elliptical galaxies is universal — 
i.e., not a function of the other observables — at least over 
the range of galaxy masses covered by the SLACS sam- 
ple. But what is the breakdown of this mass into stellar 
and dark-matter components? Fundamentally, this ques- 
tion cannot be answered without the imposition of prior 
conditions on either the form of the dark-matter density 
profile or the mass-to-light ratio of the stellar component. 
To express the ambiguity mathematically, for any stellar 
and dark-matter density profiles and pum that satisfy 
all observational lensing and dynamical constraints, we 
may transform according to 

p*^&=0p* (12) 

Pdm — > Pbm = Pdm + (1 — P)p* , 

subject only to the requirement that p^ > and p^M -> 
everywhere, without altering either the observable quan- 
tities or the internal dynamical consistency of the sys- 
tem. Previous papers in the SLACS series (Paper II, 
Paper III, Paper IV) have approached this problem from 

14 Gravitational microlensing, in contrast, can distinguish be- 
tween mass distributed among point-like stellar-mass objects and 
mass distributed in a s mooth , conti nu ous component . See , 
e.g.. IDeguchi fc Watson] lfl987l. [1981) : IWebster et al.l (H99l|): 
ISeitz. W ambsga nss. fc Schneideii U9M ) : IWitt, Mao, fc Schechter] 
( 1995): iLewis fc Irwinl l|1995l . 119961 ): ISchechter fc W ambsganss 



both the dark-matter profile and stellar-population an- 
gles. With forthcoming multi-band photometric data, we 
will be able to estimate stellar masses through more de- 
tailed stellar-population modeling, which we may in turn 
relate to lensing and dynamical masses to establish over- 
all trends within th e population, (see. e.g.JGallazzi et al.1 
2006; iBundv et al.ll2007l ; iGrillo et al.ll2008h . 

Qualitatively, our result is in good agree- 
ment with the recent th eoreti cal work of 
Bovlan-Kolchin Ma. fe Quataertl (|2005l ) and 

Robertson et al.l (|2006h . Based on the analysis of 
numerical simulations, these authors argue that the 
tilt of the FP is the result of a systematic trend in the 
central dark-matter fraction that is established and pre- 
served though the assem bly of successively more m assive 
spheroids via mergers. iRobertson et al.l (|2006l ) argue 
further that this trend is originally established through 
the importance of dissipative gas processes in disk 
galaxy formation and early merging . Similar conclusions 
were reached observationally by iPadmanabhan et al.l 
(|2004h with the aid of stellar-population modeling and 
the oretically-motivat e d dar k-matter halo models, and 
by ICappellari et al.l ([2006) through the combination 
of spatially resolved dynamical modeling with stellar 
population modeling. 

9.4. Final thoughts 

In conclusion, our results offer the following physical 
picture for the tilt and tightness of the FP for massive 
early-type galaxies (see also the discussion of Paper III) . 
The luminosity profiles are, on average, well described 
by a de Vaucouleurs model luminosity profile. This pro- 
file is embedded in an on-average scale-free isothermal 
(p oc r~ 2 ) total mass density profile, with the distribu- 
tion of the dark matter dictated by a bulge-halo "con- 
spiracy" to establish the isothermal total profile. The 
conspiracy is not far-fetched if the isothermal profile is 
in fact a dynamical attractor for the evolution of cen- 
trally condensed self-gravitating particle systems, since 
stars and cold dark matter together constitute a single 
collisionless compo nent from the point o f view of gravi- 
tational dynamics (|Loeb fc Peebles 2003). The galaxies 
thus form a two-parameter family, indexed by the effec- 
tive radius of the stellar distribution (the only observ- 
able scale length) and the central velocity dispersion of 
the system (the only parameter of the mass model) . Lu- 
minosity (and hence mass-to-light ratio) then varies sys- 
tematically with the scale size of the stellar distribution 
in such a way as to give a gradual increase in the total 
mass-to-light ratio with increasing total mass within the 
stellar effective radius. 

Observationally, this picture remains ambiguous as to 
the systematic breakdown of mass into stars and dark 
matter along the FP, a deficiency that can be partly ad- 
dressed through analysis of forthcoming multi-band HST 
imaging (ACS, WFPC2 and NICMOS) of the SLACS 
lens sample, which will give better insight into the role 
of stellar-population effects in establishing the scaling re- 
lations examined here. In addition, spatially resolved 
kinematic data are being obtained for the SLACS lenses 
dCzoske et al.l l2008h to gain the tightest possible con- 
strai nts on the possible luminous / dark decompositions 
(e.g.. iBarnabe fclCo opmans 2007}). This ongoing analy- 
sis will lead to a clearer and more robust picture of the 
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empirical physical relations described in this work. 
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